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STATUS OF STUDIES ON BUMPER PLATE PENETRATION AND SUBSEQUENT 

PHENOMENA 

SURVEY OF SCOPE OF THE PROJECTS 

Two closely associated projects are being solved with special consideration on the 
manner in which they interlock. The overall project is to obtain a computer program for 
the penetration of a thin plate of aluminum, i.e. o bumper plate, by a sphere of rock, 
i.e. a quartz sphere. In addition to the penetration of the plate, the solution is to pre- 
dict the shape of the ejected cone of plasma with entrained aluminum fragments and 
perhaps some quartz fragments. A semi empirical relation Is sought to predict the ab- 
lation of the fragments to nothing as the fragments penetrate a stack of sheets of paper. 
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Finally, the cone of plasma impacts on the stack of poper sheets and the shock pressure 
is attenuated by the counter flow of the degradation products from the decomposition 
of the paper by the hot plasma. 

Phases from both projects are reported. One project covers the initial impact, 
the crushing of the sphere of rock, the break up of the aluminum sheet and the con- 
version of the sufficiently shock -compressed regions of rock and aluminum into a plasma. 
The other project considers the ejection of a cone of plasma with entrained particles 
from the impact zone, it expansion as it traverses a region of free space and its im- 
pact on a stack of paper sheets. The ablation of fragments in penetrating a stack of 
pc ^r sheets is a part of the last project. With this interlocking of the problems, the 
key study relates to the penetration of the plote and the ejection of cone of plasma. 
Problems in other phases must await the solution to this specific phase. This specific 
phase is a part of the proposed M.S. thesis for Mr. Mcrk Hooker. 

At the initiation of these projects, a computer program was available which gave 
good results for the formation of a cro.er by the hypervelocity impact of a sphere of rock 
on a thick slab of aluminum. It was proposed to adapt the solution for the crater forma- 
tion to obtain a computer program for the first part of the penetration of the thin plate. 

It is believed that a computer solution cannot be wirtten, without very speciol corrections, 
for the entire impact which includes the penetration of the bumper plate and the egress of 
the products of the impact on the bumper plate. The reason is that the aluminum of the 
bumper plate is ruptured by tension. At the instant of rupture, the continuum foils and 
all of the equations fail since they are only valid in a continuum. After a "hand" correc- 
tion, the computer solution may proceed. 

Both projects were tentatively assumed to require two years for a complete solution. 
This final report is written midway of the study, and is, of necessity, more of a long 
progress report than of a final report. Without reference to the separation of the projects, 
the overall study may be roughly divided into the following phases. These following 
phases are identified with the phase in this report which treats this subject. 

(1) Adaptation of ’he available computer program for the formation of a crater to 
a program for the start of the penetration of a bumper plate. This phase is in progress 
and the program is discussed under Phase 3 of this report. 
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(2) The technique for making the "hand" corrections when the continuum is 
ruptured. This subject is under intensive study and a proposal is being formulated 
to continue the work next year. The discussion in Phase Three treats this subject, 
but does not enter into the subject in detail. 

(3) After sufficient compression up the Hugoniot curve, the adiabatic expansion 
will leave sufficient energy in the material so the material is properly classified as a 
plasma. The material which egresses from the impact zone is a plasma with entrained 
pieces of aluminum and rock (quatz). The expansion of a high energy plasma is dis- 
cussed in Phase Two of this report. A different equation of state is required for a 
short density range from solid density to 1/100 to 1/1000 of solid state density. After 
the insertion of this correction, which is available, the plasma calculation can be 
extended to solid state densities and to much lower energies. 

(4) The computer program for the formation oi a crater requires many corrections 
for the penetration of a thin plate. Many of tl.e corrections have been made. It is 
believed at this time, that corrections which arise from the interactions at the boundary 
between the aluminum and the penetrating rock are the most important ones which remain 
uncorrected. The work that is reported in Phase Three presents the present status of 
work on this program. 

(5) The plasma that is ejected from the bumper plate contains many fragments of 

aluminum and rock* The fragments are usual!) around the circumference of the plasma. 
One phase of a contract is to derive a semi -empirical equation to represent the penetra- 
tion of the fragments into a stock of sheets of a special paper. After a second trial at 
the derivation, an analytical equation was derived and compared with the results of one 
experiment. Since there ore three mojor constonts in the solution, the evaluation of 
all three constants is not possible with the results from one experiment. The results of 
many more experiments are available, but the author only accepted the results from one 
at the start of this program. AAany more experiments and some fitting of equations is 
necessary in order to obtain a semi -empirical equation. This work is considered in 
Phase One of this report. , 
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In the following report, there are three major phases discussed. These three 
phases do not exactly coincide with the above breakdown Into five parts that are 
discussed in this introduction to this report. In the preceding list of five parts, the 
phase in the report to which the discussion most nearly relates is mentioned. The 
phases in the report are Ablation of a Solid Sphere by Penetration of a Stack of Sheets 
of Paper, Expansion of a Sphere of Plasma with Different Initial Energies, and Adapta- 
tion of the Computer Program to a Solution for Meteoroid Bumper Penetration. 
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PHASE ONE 

ABLATION OF SOLID SPHERES BY PENETRATION OF A STACK OF SHEETS OF PAPER 

INTRODUCTION 

In the iSird Quarterly Report^ the complete ablation of a sphere of rock by penetrating 
sheels of paper was considered. The force to decelerate the sphere was assumed to have the 
form of Stoke *s law for the frictional force on a sphere falling through a viscous liquid. The 
solution was obtained, with some difficulty. On comparison with experiment, the coeffi- 
cient of viscosity was found to be ridiculously low. It was less than the viscosity o^ water. 
Stoke 's law ossumes that the losses are distributed through a substantial volume, and the 
affected volume extends fiom the surfoce of the sphere for some distance into the paper. The 
obove result was interpreted to mean that the actual losses are very close to the surface of 
the sphere. Or this assumption, a new equation was formulated for the forces that op| ose 
penetration of the sphere through the paper. 

In this report, the equo‘ion for losses near the surface of the sphere is restated with a 
little explanation of the basis on which it was formulated. A solution is eventually obtained 
in terms of three basic constants. With the results from one experiment, it is not possible 
to evaluate the basic constants. It is feasible to evaluate a ratio between two critical con- 
stants and to evaluate this ratio over the range of values which are assumed to be of in- 
terest. 

The derivation of this equation is given in detail. If there is sufficient financial support 
for the author to complete the problem, the sponsor will have no difficulty in employing 
experimental results in order to ascertain the merit of the equation. At this time, the author 
wishes to acknowledge the exceedingly helpful assistance of a graduate student, Mr. Mark 
Hooker. Mr. Hooker is working on a thesis to fulfill part of the requirements fc. his M.S. 
degree. His assistance to the author was most valuable in checking equations and signs. He 
WO:; also most helpful in discussing the study. In addition. Hooker Is entirely responsible for 
the technique by which the variables were separoled in order to calculate u In terms of r. 

The author's original concept of the plug, ahead of the sphere, was that It had the shape 
of a cylinder. As shown in Figure 1 , the plug is really a very thin shell that is ahead of the 
sohere. From thij difference, there should be a numerical constant inserted in the first term 
on the right of the equality sign. Instead of the term r,r^ /r u for the mass of a cylinder, 
the term should read 2 v r' n\ for the mass of a thin hemispherical shell . 
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STOPPING OF A SPHERE BY PENETRATING PAPER 


In the third progress report, It was shown that an absurd!/ low value is obtained for 
the viscosity, I], when the drag on a penetrating sphere Is represented by Stoke 's law. The 
low value of T] indicates that the actual losses ore concentrated very near the penetrating 
sphere. In the derivation of Stoke 's law, a particle in linear motion is restrained by 
viscous forces in o liquid with laminar, or "streamline" flow. This distributes the frictional 
losses through a considerable volume that surrounds the path of the penetrating sphere. The 
low value of T| from the application of Stoke ‘s law to the experimental data Indicates that 
the loss cannot be In a large volume but must be In a small volume and this requires that the 
losses be very close to the penetrating sphere. 

In order to confine the loss to close proximity to the sphere, a new relation is sought. 
In this relation, the losses are represented by two terms; (1) the force to accelerate the ma- 
terial in a short column ahead of the spf. re to the velocity of the sphere and (2) the force 
necessary to move the sphere and column against the restraining, frictional type forces that 
surround it and any force that is required to shear the material. The first term depends on 
the instantaneous projected area of the sphere. In contrast, the second term is the shear and 
friction of moving the accelerated material into the stack of paper. In words, the equation 
may be written 


Force to 

decelerate sphere 


Force to accelerate a 
a column of thickness 
-t and area n r'* 


Force to shear and 
plus force a plug through 
the hole 


nr^Du) 


d t 


_ d(TT r^ -tr oi u) 
dt 


2TTr ^,r S 

dt 


4 

In the first term on the left of the equality, j nr''’D is the Instantaneous mass of the pene- 
trating sphere, vher. o is the density of aluminum. When the mass is multiplied u, this 
term becomes the momentum. The differential of the momentum is the force in dynes. 
Similarly, T r' (,r is the mass of a circular plug of the paper material that is accelerated 
ahead of the entering sphere. It is believed that thi: material piles up to an undetermined 
thickness, f r, and then permits the sphere to squeeze through it with added friction which 
is given In the last term. By multiplying the first term on the right of the equality by the 
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velocity, o, the mornentum is obtained and the rate of change of momentum v/ith time 
is a force in dynes. 

The final term on the right is the effective frictional force as the circumterence 
of the sphere is forced through the plug cf material that is formed according to 
the first term on the right. The force in dynes is obtained by multiplying by the accelera- 
tion. The term, S, is the force in grams/cm‘^ of wall area of the plug which has a cir- 
cumference of 2 TT r and a length of -tr. 

The variables in the equation may be obtained in a more generally applicable term 


jP.u3r + 3 npt gy =ntp,3r „- + n«p,r ^+2"iSr 


Partially collecting terms 

u 


p - TT -t Pi j 3r‘ ^ ^ J ^ Pi j ^ ^ ^ ^ ^ 


Continue to collect terms 


j n p - TT-f, Pi ) 3r' 


, dr _ / 

^Tt i 


rr t Pi + 2 n-f, Sr*' - j np 


du 

d? 


Solve for ^ 
dt 


dr 


n -t p 1 r'"’ + 2 TT -f- S r‘" - ^ n p 


TTP * 3r‘ u - Ti-tpi3r‘'u 


j.p . 


du 

dT 


This equation may be rewritten in a form for use in the calculations. 


^ ^ + 2fS - ^ Pr ^ ^ 

dt 4 p - 3-tPi u dt 


Solution of Simultaneous Differential Equations for Relation between u and r 

With the differential equation that is derived above and with the \v/o others that 
were derived In the third progress report, there are three differential equations to solve 
simultaneously for the three Independent variables r, u and Q. These simultaneous 
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differential equations are collected here. 


■^=4TTr‘D (n 

dt 

This equation comes from the definition for the lir flow of heat. The term-^ Is 
the flow of heat in calories per second. The term, D, is an * adjustable constant'* 
which replaces the term 



where K is the thermal conductivity of the bunched and compressed material which is 
just outside the surface of the sphere. The term ^ is the effective temperature gradient 
in This material . 

The second differential equation was derived In the preceding section. This equation 
was derived from the forces acting on the sphere on the assumption of "localized forces". 

It has the form 


d£ ^ ^Pir +2tS - or ^ ^ 
dr ^ 0 - u dt 


( 2 ) 


!n this relation, r is the instantaneous radius of the sphere that is penetrating the paper, 
is the effective thickness < *f the "plug" of paper pushed ahead by the sphere, S is the 
shear force which opposes tne movement of the plug, o is the density of aluminum and 
is the density of the paper. The instantar^eous velocity at any instant is u. 

From the conservation of energy by equating the initial kinetic energy to the losses 
of energy, the following equation was obtained in units of calories per second. 


- 1.6685 X 10““ 0 (ar ^ + 2r ’ u (3) 

In this relation, the values of p,r, u and t are the same as in the preceding equations. 

Eliminate variables in order to simplify the equations. Start by equating equu- 
tions (1) and (3) in order to eliminate 
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4nr D = - 1.6685 X 10‘® o 


+ 9 '= 

+ 2r u -r 


This equation will eventually be solved for ^ and ^ . As the first step, the equation 
may be formed into the sum of two differentiols. 

o : n dr T du _ 4rr * D 

dt “ dt ■ 1.6685 k 10~"p 

Divide this equation by 3r' u and transpose one term to the right side of the equolity 

^ 4 T7 D 1 ^ L ^ /A\ 

dt 1.6685 X 10“^ p 3u^ 3 u dt ^ 

Among the three simultaneous equottons for solution, there is another differential equa- 
dr 

tion • This is Equation 2 above and is rewritten here for ready referenc - 

dt 4 d-'3^x udT 

The terms on the right hand side in Equations 4 and 2 may be equated to each other 

since both are equal to ^ . 

dt 

t r + 2 AS ^ I P r _ 4nD 1 2rdu 

u dt 1 .6685 X ” p 3u‘ 3 u dt 

Rearrange the terms in order to factor out ^ 


A o^r 2^S - ^or 2 \ ^ - 1 

\ 4 0 - 3£pi 3 ' u dt 57^ 1 . 


4 TT D 

6685 X 10“ ” 0 


In order to simplify this relation, multiply through by -3. 


(4p - 3Ap^ ) r - 6 A S 


In - 3 f D 


6 A S - \ du 
; ^')“dT = 


4 TT D 

1.6685 k I0“*p 


In order to further simplify the relation, introduce the constants Pj, Pg, and P;, where 
the constants are defined by the following relations 


1.6685 X 10 


Introduce the new constants and further simplify 
/ rP= + Pn du _ p. 


trj5;Pp ~4p - 3A Oj; = - 6^S 


-2r u^=P, 
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Obtain a common denominator and rearrange the terms 


( 5 ) 


du . 


This is an equation for in terms of the variables u artd r and of the constants , Ps 
and P ^ . 

The next step is to use the same two Equations 4 and 2 in order to obtain an expression 
for P/oceed b/ obtaining an equation for from Equation 4. 


2 ^ ^ ^ _ 1 p 

3 o dt dt "37*’ ' 


du _ 3u dr 3u 1 

dt 2r dt 2r 3u“ 


(a1 


Now consider Equation 2 and modify it in order to obtain an equation with the constants 
inserted. With a slight rearrangement. Equation 2 becomes 

^ ^Oir ~^pr-t2^S ^ _ l_ rPg + P3 ]_ ^ 


dt -3 


4o - 34oi u dt 


u dt 


Rearrange the above relation In order to solve for ^ 

dt 


^ = -3u_Lb_ * 

dt ^ rPa+Pa 

In order to eliminate equate Equation a to Equation b. 
dt 


3 u d r 3u 1 


P, =-3u 


Pa 


2rdt 2r 3u" ■ ‘ rPa + Pj, dt 

Divide both sides of the equation by -3 and rearrange terms 

/, 1 Pq \dr 1 P 

\ 2r ‘rPa + Paj dt ~ ~ STT’ 

Continue to simplify this relation 


(b) 


-3u® * =P, 

P,+Po dt 


( 6 ) 
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Starting with Equations 4 and 2, ^ is obtained in Equation 6 as a function of the variables 

du 

u and r and of the constants Pp and P^. In Equation 5, — was obtained as a function 
of the same variables u and r and of the same constants. Equations 5 and 6 were obtained 
from Equations 4 and 2. No variables or constants have been eliminated but the equations 
arc in a different form. 

From an inspection of the relations in Equotions 5 and b, each equation may be solved 
for the constant, P^ . The two equations may then be equated to each other by the elimina- 
tion of the variable, P^ . 


P^ - rPp du _ o : Pn ~ 

Pp *^dt Pr,+rP:= dt 


After simplifying this relation, the following form is found 

1 du _ 1 dr 

3u dt Pn r dt 

p ^ 

' p 

In this relation, the variables are seporoted and integration may proceed at once 



1 

1 


u 


r 




P p 


A relation between u and r is necessary in order to evaluate measurable quantities in the 
experiment such as the length of path to stop the penetrating particle, or the length of 
path to completely ablate the sphere. 

In the above relation between u and r, a constant of integration should be added. 
From the manner in which the equation was integrated, the constant could be inserted in 
several ways. In one method, the equation may be rewritten as 


1 

^Cu)^ 




Li 

Pp 


In this case, the constant, C, may be considered as a multiple of the velocity, u. There 
will be no other reference to C until the discussion of the evaluation of the constants. 
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LENGTH OF PATH THROUGH THE PAPER TO STOP A SPHERE 
(Duration of Period of Penetration) 


From the relation between u and r that was obtained in the preceding section^ the 
length of path through the stacked shot ts of paper may be calculated. There are two solu- 
tions for the length of path: (a) The length of path for complete ablation of the entering 
sphere. This occurs when there is an excess of the initial kinetic energy, (b) The length 
of path through the stacked sheets of insulation to stop the sphere with a finite, final radius. 
No entering sphere in the experimental range of sizes is expected to reach the rear 
surface of the stack of sheets of paper. If the solution from this study proves to be adapt- 
oble to obroin reliable results, then it would be desirable to obtain a solution for the 
velocity and rodius of a sphere which is either oversize and has excessive initial energy in 
order K predict the radius of the sphere and the velocity at the back plate after It has pene- 
trated the insulation. 

In order to make the desired computations, it is desirable to select three equations from 
the preceding section as the foundation for the studies in this section. The first of the three 
relotions is the relation for ^ when ^ is eliminated from the equation 


P3-rP, 

Pp 


du _ 

dT ' 


Px 


(5) 


The second of the three equotions is the relation for ^ when ^ is eliminated from the 
equation 


^ : P.T -rP, dr 

f-rrf df 


( 6 ) 


The final equation of the three equations Is the relation between r and u which is obtained 
from the two preceding equations 


1 


u 


I 

n 


r -t 


Li 

P 2 


(7) 
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Length of Path of Sphere in Stoclced Sheets of Paper 

From the three Equotions 5, 6 and 1 , it Is desired to calculate the length of path of 
the penetrating sphere In the stacked paper insulation. There are two distinct problems for 
wb'ch a solution Is required. The length of path is desired (a) to stop the sphere by ablating 
it o nothing and (b) to stop the sphere by ablating it down to an approximate sphere of finite 
rv.dius, r^. Start with Equation 7 but rewrite It in a form that gives the value of u 



ds 


Th s relation may be written in another form by recalling that ^ where A s Is an element 

of length along the path of penetration of the sphere. The new form of the above equation is 


\ / dt 


= 1 


Introduce a new independent variable, r. 


/ P- \'* ds dr , 

rpt) d7^-' 

This equatian may be rewritten in still another form 


Pa'" 

\ Pa 

dr . 


, \'’ dr _ dr 
) dt ds 


The value of is known from Equation 6 and on substitution, the relation becomes 


' 4Pn\7 h\(LLlilA =± 

; ^ p; ) rJu^jVp.-, +rpj ds 
Now insert »’ie value of u‘ from Equation 7 




I ' 3' • ■ ] ( -la j( PaLiZaV + P.-' ) 

I Pa A 3 A I (Pn-rP.) 


GImpl'fy the above relation by collecting terms 


p, 


dT 


(P.-r -rPa) \ 3P^'[ 

separate variables and put in a form for integration 


ds 
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(8) 


Since the preceding integral is rather Involved, the terms on the right are inte- 
grated one by one. After integration, the separate solutions are reinserted into a single 
equation. Consider only the terms in the fractions that are involved in the first term on the 
right in the last equation above. The expression to be integrated is 

dr 

From a table of integrals, there is a substitution for this type of term. The substitution 
requires that y " o + bx, then tables of integrois show that the following relation exists 




For the integral under consideration, m = "10, a b = Pg . With these substitutions 



(P3 + Pprr^^d(P3 ^P^r) 


1 1 
9P^ (P^Tp;^ 


The involved terms in the second expression for integration ore 



From a table of integrals, there is a substitution when the following identity hold for this 
type of equation, y ~ a + bx. The relation from the table of integrals is 


/ x^ dx 
(a + bx)''^ 



For the integral under consideration, the following relations hold: m = 10, n= 1, a = Pj, 
and b = Pg 

With these substitutions 
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J -8 

« y 


8Pa 


9 ? 


-3- y 


-9 


_ 1 1 

■ (Pa + Psrr 9P71Prrp7p 

Now coliect all of the terms and reinsert them in Equation 8. For an integration without 
limits, such as this one, a constant of integration, C, must be added after the simplifica- 
tion is completed* 


3 ^ 


s 


Pn 1 1 * Pi > 

9Pa (Pa + PsO" ^8Pa (P. +P.r)" '9P, (P, + P^r)^ 


Simplify by several steps, change signs through the equation, multiply by Pg, and com- 
bine terms 


Pi - Pa n . n 1 1 

W* * (Pa + Pa r)» [9 9j ' 8 (^ + PaO^ 

Put the equation in its final form and add the constant of integration, C 


Pi _2P.a ^ 1 

9 (P~iT7P "8 (Pa +P 2 r)“ 


( 9 ) 


The preceding relation. Equation 9, is a general equation with the length of 
path, s, and the radius, r, of the sphere* This equation may be evaluated between 
several limits. Two sets of these limits are of particular merit In evaluating the unde- 
termined "constants" that are in this equation. One set of limits corresponds to a sphere 

of radius, r , that enters the stack of paper at the start of the path of penetration, s = 0, 
o 

and is completely ablated, r = 0, after a path length, s^. The second set of limits corres- 
pond to a sphere that is not completely ablated and stops in the stack of paper. This set 
of limits corresponds to a sphere of radius, r^, that enters the stack of paper, s = 0, and 
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is ablated to a radius, d , after traveling olong the path of penetration for a distance, 

s^. In other words, a portion of the entering particle is found in the stack of paper at 

a distance, s ,, from the entrance of the particle into the stack of paper, 
d 

The first set of limits, which are to be inserted. Is for the spherical particle which 
is completely ablated os it penetrates the stock of sheets of paper. The limits for inte- 
gration for this set of conditions are: initial conditions, s = 0 when r ” r^; and the final 
set of conditions, s = s^ when r = 0. When limits of integration are inserted, there is 
no need for the constant, C. Omit this constant and insert the above limits of integration 
Into Equation 9. 




S 


S = s. 


s = o 


2P, 1 

9 (Pn +Psrf 


r = o 


r ~ r 


1 1 

8 (P, +P,r)“ 


On the insertion of these limits, the following relation is obtained. 


= o 


r 

o 


Pi „ s =21^ r 1 
3(prr f 9 [(p^ 


1 

(P r» ^ P 



1 


(P.+ P 


ro)«. 


Separate the terms with r^ from the other terms 


1 Pn + P-Tq 2P^ 1 ^ 2Pg 11 

8 (P.-, + Pi ro>' ‘ 9 (Pa + Paro)^ ~ 8 (P3) 


Collect and simplify the terms In this relation 


(Pa+Pi^ro)" [■ 72 " 8 ^ (P;? 


Solve for Pj s^ and continue to simplify the terms 

P . = (Lil 9P: ro-7P.., 7 _(^ 

f 24 (Pa +P;jro)^ 24 V p J 


(10) 


To recapitulate, this solution is for a sphere of radius, r^, which enters a stack of paper 
sheets at a length of path through the paper of s =0. After traveling through the sheets 
of paper for a distance, s^, the sphere is completely ablated. For comparison with this 
solution, derive the other practical example in the next paragroph. 
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Consider the other set of practical limits when an entering sphere is stopped in 
the stack of paper before It Is completely ablated. The limits of integration for this 
set of conditions are 


s = 0 when r = r 


s = 5 , when r = d 
d 

In this set of relations, s , Is the length of path through a stack of sheets of paper before 

d 

the sphere is stopped with a radius, d. The integral with the limits of integration becomes 

r = d 


r - r 



s - s , 

r = d 

A s 

1 

1 1 

3(P2V 

9 

s — o 

^ ^ 8 (P,+P2r)® 

o 


With the limits inserted, the relation becomes 


*d’“r 


1 1 
(P,+Ppd)^ ■ (P.,+ Par 




+ Prd)“ “(Pa+P 




Separate these terms into two groups according to d and r^, which are the final and the 
initial radii of the penetrating sphere, respectively. 


(P^ 




2 ?^ 1 


- - P - — P 


Continue to simplify and collect the terms in the above relation 


_ TP, - 9Pgd 1 7P^ - ^Ppro 
72 (P^:, + Pa5^ ■ 72 (P>Par')^ 

Solve this relation for Pi*j simplify the terms 


n . -(Pz)° TP, -9P,d (Pg)® 7P,-9Pa 

^'*d 24 (Pa + Ppd)'' 24 + P'proT 


( 11 ) 


The preceding two Equations 10 and 11 are for the complete and the partial ablation, 
respectively, of a high velocity sphere which enters a stack of sheets of paper. One set 
of experimental data is available to evaluate the constants in Equation 10 but no data is 
available to evaluate the constants in Equation 11. 
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LIMITED EVALUATION OF CONSTANTS BY COMPARISON 
OF DERIVED EQUATION WITH AN EXPERIMENT 

After the analytical form of the above equation is obtained, it is desirable to evaluate 
the constants. At this stage, over optimism must be avoided. The general equation v^lll 
probably fit a very restricted range of the data and the equation must receive a major 
modificotion befor*^ it may be applied to the general problem. Since there are three 
basic constants, Pj, Ps and P^/ these cannot possibly be determined by the results from a 
single test. (In order to keep close contact with the Sponsor, the author has only accepted 
the results from a single test up to this time.) With a single experiment, only a range of 
ratios may be obtained between the most sensitive "constants". In this section, this ratio 
is obtained for adjustments In two sensitive "constants". The first step will be to restate 
the results of the one experimental test which has been accepted by the author. 

In the third Quarterly Progress Report, the values are given for one experimental test. 
This information is for Experimental Shot No. 203. The incident particle was an aluminum 
sphere with a weight of 0.010 grams. The initial diameter of the sphere was 0.075 inches, 
or 0.1905 cm. For use in the derived equation, the initial radius, r^, was 0.09525 cm. 

The initial velocity was 6.15 km/sec, or 6.15 x 10^ cm/sec. The length of path through 
the stack of paper sheets in order to ablate the sphere was about 8 cm. The mean density 
of the paper was 0.4609 gms/cm'^ which is less than half of the density of water. 

The derived equation employs three constants, Pj , Pg, and P^ , These constants are 
defined in the preceding derivations. There are actually more than three constants, but 
for the following evaluation, it will be found that the most significant adjustable constants 
are in Pj and in P 3 + Ps ro - It is obvious that a single measurement cannot evaluate 
these three constants. It can only approximately determine a ratio between the two most 
sensitive "constants" among those which determine Pj, P^ and P 3 . 

The two most sensitive "constants" are and 0. The active constant, £S, is the 
product of times a combination of the average shear and the average friction. The method 
of combination of these two factors Is not known. The value, Is defined by means of the 
product where this value Is the effective length, parollel to the direction of mot on, of 
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the friction surface around the circumference of the penetrating sphere. In the product, 
lx, <’ is a pure number between 0 and 1 and r is the Instantaneous radius of the penetra- 
ting sphere. In other words, ^S, is a combination that is proportional to the fr'ctlon and 
to the shear in unknown proportions. The sensitive "constant", D, Is representative of 
entirely different factors. It was discussed extensively in the preceding Progress Report. 
It is defined by the equation 


D 



where K is the thermal conductivity and-g^ is the instantaneous radial temperature gradient. 
The definition applies to the region of damaged material which is between the undisturbed 
sheets of paper and the surface (probably molten) of the penetrating sphere. 

The equation for evaluation is Equation 10 v/hich has the form 



The constants In this equation have the following values. 


n 4 TT D , dT 

As stated above, K is the thermal conductivity of the crushed and torn material that 
surrounds the penetroting sphere. The temperature gradient, is through this same 
crushed material in the region which is just outside of the metal of the penetrating sphere. 
The numerical value in the denominator is for the purpose of converting calories and thermal 
conductivities into units which are compatible with the c.g.s. units in the other differen- 
tial equations. The manner of obtaining this numerical value is discussed in detail In the 
third Progress Report, The value of D is not known and may only be estimated. A range 
of estimated ratios will be obtained from the experimental test. 

The constant, Pp, Is given by the following relation 

Pg =4o 

The density, p , Is the density of aluminum which is assumed to be 2.765 grams per cm'^ 
from various references. The lensity, Px , Is the mean density of the stack of sheets of 
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paper. From the experiment that is referenced, the volue of is 0.4609 grams per cm^. 
This density is less than half of the density of water. 

The constant, , is defined by the following relation 
Pn = S 

There is no guidance known to the author from which to select numerical values for I and 
S. For the calculations that follow, a value of 10 grams per cm' Is assumed for S. Then 
values of i ore assumed which give values of D which appear to be In a practical range. 

The derivations of the three simultaneous, differentiol equations and the search for 
their solution has been In much more detail than is usual. This detail was added because of 
the great uncertainty as to the basis for the derivation. It will make it much easier to 
modify as comparison with experiments proceed. The insertion of numerical values is Ulus* 
trated by one solution which is reproduced in Appendix A. Three solutions have been ob- 
tained for different values of i which determine three values for D. These results are re- 
ported III Table I. 


Table I 

Values >f D Calculated from Assumed Values of 
I D 

Pure Number calories 

cm^ sec 
0.01 12.90 

0.011 103.9 

0.012 1516. 

It is to be recalled that there is a minimum of three constants; P^, Pp and P-,. It is 
obvious that a single experiment cannot uniquely determine these three constants. Actually, 
there are more than three constants and the author has chosen two of these constants to be 
more important and has assumed that they vary more rapidly than the others in range of 
particle size and range of penetrations tha;* are represented by the one experiment, for 
which data is Immediately available. The choice Is arblirc^'y and a better choice might 
be made. The entire subject is so uncertain that It is discussed in more detail in the next 
section. 
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DISCUSSION OF EVALUATION AND SIGNIFICANCE OF CONSTANTS IN THE 

DERIVED EQUATIONS 

In fhls secK.n, /andom comments are given on the evaluation of the "constants” In the 
derived equations by means of comparison with experimental data. The derived equation 
must be treated as :>uspected of not satisfying the experimental data. For this reason, the 
experimenlai data should be inserted In a way that will emphasize the deviation and permit 
the basic, differential equations to be modified in order to be more correct. As an example, 

K must be suspected of being a function of the velocity, u, or of the radius, r, v^+iich is 
essentially the same thing, In this stage of the development. Since these r nts will 
range over a wide range of subjects, subheadings are employed although he sub- 

headings may only be a paragraph in length. 

Constants Involved In the Evaluation 

For simplicity In integrating the differential equations, three constants are employed. 
These constants ore Pj, Pp and P.:, . As indicated In the preceding section, these three con- 
stants represent at least 6 constants: K, P, pj, i and S. In practical examples, P 
and Pj, are known in the most general cases without conducting an experiment f j evaluate 
them. The density of the \ snetrating sphere is P, and the average density of the stack of 

sheets of paper is . Two of the constants are impossible to evaluate separately by the type 

dT 

of measurements that are presented in the illustration in this report. These are where 

K is the thermal conductivity and ^ is the thermal gradient in a region of thickness, A r, 
which is just outside of the penetrating sphere. Minor modifications may be made in the 
experiments which will permit separation of these constants. 

The constant, i, is difficult to define. The explanation is perhaps easiest to Introduce 
by means of a sketch which is designated as Figure 1 . As the sphere penetrates, it Is assumed 
to collect and accelerate a thickness, iv, of the debris from the degradation of the paper. 

The thickness of the added materials Is ^ r where ^ Is a pure number and r Is the Instantaneous 
radius of the sphere. This thickness only covers the entering hemisphere where it is accelera- 
ted by the moving sphere. From the numerical statement of this assumption, It Is to be ob- 
served that the rhlckne'S of the layer varies as the instantaneous radius of the sphere. In 
addition, there Is a roughly cylindrical shape around the circumference of the sphere where 
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fhe crushed material that is ahead of the sphere “slides" betvyeen the largest diameter of 
the "instantaneous** size of the sphere. This is also the position at which the penetrating 
moss shears the sheets of paper. This cylindrical shape is a friction surface as the sphere 
penetrotes. The combined shear and friction is tentatively assumed to be expressed as S 
grams per cm*. An entirely arbitrary, numerical voiue is assumed for S, since It is necessary 
on account of the manner In which the Initial evaluation of the constants is mode. The value 
of S is orbitrorily assumed to be 10 grams per cm . This value may appear to be small but 
recall thot the overage density of the stack of paper is 0.4609 grams per cm which is to 
be compared to 1 .0 for the density of water. In addition, the surface of the penetrating 
sphere is at the mi nmum temperature which corresporKJs to the melting temperature of 
aluminum. 

Omission of up and u^ from the Evaluation of the Constants 

In the numerical evaluation of the constants in Equation 10 by the Insertion of the results 
of an experiment, all of the experimental data is not employed. The values that are reported 
from the experiment are the Initial radius, r^, of the penetrating sphere, and the final rodius, 
0. The density of the penetrating sphere is o and the average density of the stock of poper 
is . The length of the pa.S, s^ for complete ablation is also employed. Neither Equa- 
tions 10 nor 11 contain u, explicitly. This is a consequence of the relation that was de- 
veloped between r and u. This relation was obtained from the simultaneous, differential 
equation from which the solution was eventually determined. 

To obtain the final differential equation for Integration, the velocity, u, was eliminoted 
from the differential equations for integration. The velocity, u, was eliminated by using 
the relation that 'was found between u and r. As a consequence, the velocity, uq, does not 
appear explicitly in either Equotion 10 or 11 . From a single experiment, it is olmost im- 
possible to tell much about the velocity, u. By reference to the poge where the relation 
between u and r is obtained, it is seen that the most general solution requires a constant, C, 
which may multiply u. The constant could be inserted in another place in the equation, 
but this position is perfectly general. Until the constant is evaluated, the substitution of 
u for r in Equation 10 will only evaluate the jstant, C. This constant may, of course, 
be evaluated more easily by a direct substitution in the Serw. en u^ ond r^. The 
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other apparent limit, the final radius, r = 0, and the apparent final velocity, =0, 
connot be employed for the evaluation of C. The reoson is that the following relation 



is indeterminant when u is zero. Division by zero is forbidden. It is very logical that 
u should rK>t be employed when r = 0. It is not necessary that the velocity should 
approach zero at the same time thot the radius of the sphere becomes zero. 

By reference to the equation between r arKi u, it is interesting to use the values for the 
constants, arui with the ratios that are reported in Table I. With this substitution, 
the opproximote values for C are given in Table II. The equation that relates u to r is 
given just above this paragraph. The following. Table II extends the results in Table I 
in the preceding section. 


Table II 

Values of C and D Calculated from Assumed Values of 

C 

1 

D 

Pure Number 

Pure Number 

calories 

1.112 

O.OI 

cm' sec 
12.90 

.727 

0.011 

103.9 

.441 

0.012 

I6M. 


Evaluation of Constants with Restraints on the Variables 

The next step in the study of the equations for the depth of penetration is to evaluate 
the constants in a way which will yield the most information on the nature of any possibly 
omitted factor. From this information, the initial, three differential equations may be modi- 
fied, or another equation may be added. The practical evaluation of the constants will 
probably depend on the experiments which have already been made. Since the number of 
experiments for the evaluation of the constants is probably limited, the evaluation must 
proceed in a flexible manner while keeping the general objective in mind. The objective 
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Is to evaluate groups of measurements In different ranges of the variables. The equation 
Is good provided the short range particles give the same values for the constants as the 
tong range particles. 

To a limited extent, the variables Tq, Uo and may be considered as the most probable 
variables wf»ich influence the secondary constants. The secondary, very critical constants 
are represented by D, t rnd S. If this selection is found to be correct, try to vary the 
voriobles over a limited range in order to emphasize the deviations in the values of the 
constants. As an exomple, consider one of the principle variables to be the length of 
path, s^. The observed lengths of path could be divided Into four ranges; such as short paths, 
short intermediate paths, long intermediate paths and long paths. Use all values of Uq which 
give path lengths in one selected range. Evaluate the constants from the data for each range 
of Sp If the equations were correct, the constonts should be the same for each range of Sp 
Then consider the same groups for the lengths of path, Sp but select the intitlal rodii, r^ , 
of the penetrating sphere. The comparison of the difference in the results for ond r^ 
should indicote whether the relation bet^en r and u is acceptable. 

It should be emphasized that this manner of evaluating the constonrs is not the best way 
to obtoin constants for the entire ronge of application of the basic equations, provided that the 
equations have been established to be correct. 

Discussion of the Values of the Constants as Determined by One Experiment 

There e three, grouped constants for determination, Pj, Pp, and P 3 . These three constants 
canfx>t be evaluated from a single experiment, because there is not sufficient data to evaluate 
three constants. As indicated elsev^ere, three representative ratios were determined between 
the Important secondary constants, i and D. The evaluated ratios are in the range that the 
final values for these two constants is expected to occur. The values of ^ are In the range from 
0.01 to 0.012. These values indicate that the thickness of the collected layer of degradated 
paper, ahead of the sphere, varies in thickness from 1 percent of ^*^e radius of the sphere. 

From the differential equation that was derived, starting with the first page on this phase, the 
first term on the right side In this equation is the rorce to accelerate this layer of paper frorr 
at rest up to the velocity of the sphere. This term appeors to be relatively negligible. The 
really significant part of is in combination with the constant S, which represents the com- 
bined friction and shear in the second term of the differential equation. 
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The really critical nature of the value of f appears in the term, (P + Parp). 

By reference to the definition of the conslants, P^ - -6 f S and Pg = 4;) - 3^Pi. 
From the example for the evaluation of the constants that is given in Appendix A, the 
numerical values for evaluation are 


7.5313 X lO^Ds 


_ (P.)® 

f ■ (p; + 


0.56955 + Negligible 
term 


(p^)® 

"(-6/S + 1.052)' 


0.56955 


When S = 10 granrj per cm^ and i - 0.01 , the term on the right becomes 


(11,062)'^ 
(-0.6 + 1.052)^ 


0.56955 


From the numerical values, it is apparent that the length of path s^, in the stack of paper 
varies as cn inverse function of /. 

There are three coses for which the relative values of -6 £ S and Ppto should be con- 
sidered. In the first case, the value of -64 S is very much less than P^ro . In this particu- 
lar case, the penetration of the particle into the paper will be relatively small. In the 
secorKl case, the term, -6 4 S, is initially smaller than Pgro but approaches the value of 
Pgro . With this range of values, the depth of penetration increases at the start and then 
increases more and more rapidly toward infinity as the two values approach closer and 
closer. Finally, the third cose, the value of -6 4 S > Pgr^^ . In all of the examples in 
this case, the length of path, s^, is negative and ihis case does not represent a legitimate 
application for this equation. 
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PHASE TWO 

EXPANSION OF A SPHERE OF PLASMA WITH DIFFERENT INITIAL ENERGIES 

INTRODUCTION 

In the introduction to this Final Report, it was indicated that progress was to be 
reported on four phases. Two of the four phases of the work are concerned with the 
problem of the penetrotion of a bumper plate of aluminum by o spherical stone of quartz. 

This phase. Phase Two, is concerned with the expansion of a sphere of plasma. The 
third ond next phase of this report will discuss the progress on the early stages of obtain- 
ing o computer program for the penetration of a sphere of rock through a thin plate of 
aluminum. As the penetration by the rock proceeds, a shock front compresses the volumes 
of aluminum and of rock in the curvilinear squares by different amounts. Every shock 
compression increases the entropy of the compressed material. In the most general cose, 
the expansion after compression is along an isentrope, i.e. along o line of constant 
entropy. This general case would require some energy to be dissipated during the expan- 
sion. Metals thar are shock compressed, usually are considered to expand along on 
adiobat with no dissipation of energy. This means thot after adiabatic expansion, from 
shock compression, the metal retains a considerable amount of internal energy in the 
form of heat. This brief statement is an oversimplification of several papers, such os 
the poper by G. E. Duvall and G. R. Fowles,^ and of other publications that are cited 
in this particular paper. 

The heat that is retained by the metal after shock corrpression and adiabatic expan- 
sion to atmospheric pressure may range from simple heating of the metal, to producing 
molten metal, to producing a plasma. The total amount of heat is dependent on the magni- 
tude of the initial compression. The thesis by Hcrdage estimates that shock compression 
to 37-1/2 kilobars w:!l result in leaving aluminum molten after it has expanded, adiaboticolly 
to atmospheric pressure. Shock compression into the megabar range converts the aluminum 
metal into an aluminum plasma. In the penetration of the bumper plate, the energy con- 
tent is divided into kinetic energy and into internal energy for each of the mesh volumes. 

As is customary with computer solutions, the solution is for finite differences. This re- 
quires all space to be divided into cubes, or regularly shaped volumes for cartesian coordi- 
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nates. When sphencai coordinates ore employed, as in this problem, the volumes for con- 
sideration are curvilinear cubes. For each of these cubes, the kinetic and the internal 
energy are given in the computer print-out. Since the expansion is adiabotic, the 
Internal energy of the curvilinear cube should not change during the expansion. This 
internal energy may be employed at any time tc identify the end product as a hot 
metal, a molten metal, or os a metallic plasma. As a consequence, a shift from the 
Mie-Grunesen equation of state to the plasma equation of state may be made at any 
time during the expansion. 

After a one by one change of the curvilinear cubes from compressed metal to com- 
pressed plasma, the expelled material from the impact zone must continue to expend. 

More specifically, the material that is expelled by the rock penetrating the thin plate 
must expand as a mixture of aluminum plasma, rock plasma, molten aluminum, molten 
rock, pieces of hot aluminum metal and perhaps pieces of hot rock. The plasma portion 
of the expelled material will expand to interlace the entire mixture. 

For this particular phase of the problem for study, the preceding discussion does 
not need to be followed farther from a general viewpoint. If an aluminum plasma Is pro- 
duced, it must be shown that this group has the copability to consider the exponslon from 
a dense picsma to c very tenuous plasma. The initial work on the expansion of a plasma 
by a member of thi: group was by R. E. Bruce. The work on the expansion of a sphere 
of plasma was the subject of the Pli.D, thesis by Dr. Bruce. This ex -student from this 
group is presently an Associa e Professor of Physics at the University of Texas at El Paso. 
The thesis by Dr. Bruce has been distributed to the mailing list on this contract, as well 
as is presently recalled by the author. For this reason, reference is frequently made to 
bis thesis, instead of to the original references. There are several ports of Dr. Bruce *s 
thesis which were in error. In particular, there are serious errors In the curves In the 
original thesis and in the manner in v/hich the solution behaved near to the vertical 
axis along a diameter of the sp^.ere of rock. Dr. Bruce was paid, as a consultant, to 
come to UAH ^or the purpose of teaching a present student to run his program end to re- 
run the programs in the original thesis. In order for the results to be more reliable, the 
curves in the remalnoer of this section are fiom the new runs by Dr, R. E. Bruce. 
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GENERAL BASIS FOR A CALCULATION OF THE EXPANSION OF A SPHERE OF 

PLASMA 

For an understanding of the method of calculation, reference is made to Chapter V 
in Bruce’s thesis, (pp. 77- 89) There is nothing particularly unique about the formula- 
tion of the equations. Scaler pressures are employed instead of tensor pressures. As is 
usual for a scaler pressure, an equation of state was obtained. This equation of state is 
unique as far as references were concerned at the time that the thesis was written. At 
the first of this calendar year, Bruce informed me that a similar equation of state with 
this elaboration has not yet been published in the literature that he had investigated. 

The evaluation of the accompanying curves in this phase of this report requires an under- 
standing of the equation of state so this subject is discussed in the next section. Before 
proceeding with a discussion of the equation of state, two other subjects require some 
discussion. One subject is concerned with the effect of microfields and the other is the 
equation of state at extremely high densities. 

Microfields in Dense Plasmas 

A discussion of the microfields in the plasma is essential for understanding the 
method of procedure. The subject of microfields will recur in this report, so a quick 
introduction is given here. In the plasma, there are ions of various specie which are 
intermixed with sufficient electrons to give a neutral plasma. As an example, a par- 
ticularly hot plasma may have Al^‘ , Al^, Al ions in a proportion that may be 
predicted from their relative ionization potentials in the plasma and from the applica- 
tion of Boltzmann's law. These Ions and the matching electron cloud for neutrality 
produce microfields and these mioioflelds lower the ionization potential for the pro - 
duction of ions. It is immediately apparent, that increases ir, the density of the ions 
and their accompanying electrons will increase the magnitude of the microfields and 
this further lowers the ionization potential to produce the ions. 

The lowering of the Ionization potential by the microfields will be considered 
again under the discussion of the formulation of the equation of state for the plasma. 
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The calculation of the number of each specie of Ions would require a total of 16 
simultaneous equations Instead of five and probably from five to nine times the 
computation costs. In order to avoid this addltionc I cost, an average Ionization 
is calculated. This is a peculiar type of average that will give roughly the same 
characteristics in the plasma as the true distribution of ions would give. This 
average number of ions is the ora i note for one or rhe curves !n most of the graphs 
that represent the expansion of the plasma. 

Equation of State Employed for Subsequent Calculations 

In many studies of plasmas^ the equation of state is taken as the perfect gas 
law. The perfect gas law considers the atoms to exert no forces on each other, 
except by collisions. In other words, the atoms behave roughly like a cloud of ping 
pong balls; i.e. there is no interaction except for completely elastic forces which are 
active only when the bolls collide with each other. For a complete description of 
all of the terms in the equation of state, the reader is referred to Bruce’s thesis. For 
the present discussion, the equation is presented and the nature of each term Is Identi- 
fied in very general terms. 

In the presented equation of state, there are terms which affect the plasma by 
different amounts in different circumstances, and some terms may often be omitted. 

The general form of the equation is given in terms of the internal energy which is ex- 
pressed in terms of the temperature and the density; i.e. the equation of state has the 
form - E(o, T). In every case, the plasma Is considered to be formed of alumi- 

num atoms, ions and electrons. It is to be recalled that the plasma is always considered 
to be neutral; i.e. the plasma has no ne^, average charge. There is a minor exception 
to this general statement, there may be local charges that build up in oscillations be- 
tween the positive Ions and an oscillating electron clouu. These oscillations consti- 
tute one term in the equation of state. After this preliminary Introduction, the general 
form of the equation of state is given by the following relation: 

E -E +E tE tF+E tE t-E 
‘'tot id ion EXC in DEG RAD OSC 
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In this relation, ^jqj total energy in the plasma. Each term will be discussed 

in some detail, but only in sufficient detail to Identify it. There is not sufficient detail 
to perform the calculation for each term. 

The first component of the energy is E ^ where the subscript is an abbreviation for 
ideal. This term represents the kinetic energy In the perfect gas law. This energy is given 
by the relation 

where represents the number of electrons and C. represents the number of ions of each 
specie, i. The ions include C which is the neutral aluminum atom. It is to be observed 
that each electron and each atom has the same average energy. This relation is immediately 
recognized as 


2 


m v 



The second te^m is deisgnated as This is the energy that is required to ionize 

the different specie to form the plasma and this energy must be lowered to compensate for 
the effect of the microfields. Formally, the energy may be indicated in the following 
manner 

+Cn(l,*M+ 

There is difficulty in evaluating this relation when it is expressed in this particular form. 
The average ionization is calculated by a redistribution of the energy between all of the 
different available forms energy that are expressed in the equation of state. This re- 
quires reiterative calculations for the solution to converge on the proper number of ions 
with the necessary number of electrons to obtain a neutral plasma. Designate the energy 
that Is found by this calculation to be I. Then the quantity, Is given by 


’ION 


A! 


where A 1 is the Ecker and Kroll correction for the effect of th3 micro^Ields on lowering 
the ionization potential. The technique for obtaining the average Ionization and the 
lowering of the ionization potential by the Ecker and Kroll method H outlined in Bruce's 



32 


thesis. It should be mentioned thc^ Bruce says that he finds it necessary to recalculate 
the convergence on every third cycle through the program. The last equation above Is 
given in order to emphasize the nature of the solution that is obtained and is not a 
description of the mathematicol steps that are employed. 

The third term on the right is and this term represents the excitation energy 

that Is in the equation of state. Every available energy will \r^.e electrons In it accord- 
ing to the maximum number of electrons in the plasma and the possible levels in .ne 
Ionized atoms and ‘he degree to which the excitnd levels are filled. The excitation 
levels are filled according to their energy by the Boltzmann relation. A general, 
analytical representation is 



The term on the extieme right is the Boltzman relation for filling the q excitation levels 
in the ions. There is a difficulty with this relation as solid density is approached, or as 
*^he density of ionization becomes very high. The microfields, in either case, exert suffi- 
cient forces on the ions to require quantization of the energy levels. This requires the 
application of the Pauli exclusion principle to limit rhe possible number of electrons in 
the excited states. This quantization is important and is mentioned later in this report. 

The fourth term is designated as This is an evaluation of the energy that arises 

from the interaction of the microfields between the different particles in the plasma. The 
specific name for the method of calculation is the configuration integral. 

In more familiar terminology, it is referrtd to as the calculation of the cluster Inte- 
grals for all particles, pair-wise. This is often mentioned in the literature as the two- 
body cluster integrals. The interaction energy is calculated for all particles, considered 
two at a time. The correction to the equation of state is particularly Important In dense 
plasmas but trailed off to a negligible correction at densities from 10“ ' to 10*"^ times the 
density of solid aluminum. It decreases wi^h the density until the correction becomes neg- 
ligible at the above limiting densities and the equation of state converges on the predic- 
tion of the perfect gas law. 
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The fifth term for the correction to the equation of state is where the 

DhG 

subscript refers to the degeneracy. This correction comes from the work of Wigner and 
is described by Bruce. The term varies as the order In the plasma decreases from the 
density of a solid to the density of a rarified gas. 

The sixth term is radiation energy that is generated and is probably 

not trapped in the hot plasma. For plasmas in the temperature and density range which 
are of interest for this project, this term is neglected. By neglecting the energy in the 
radiation, the energy in the sphere remains in the sphere. If the energy v/ere allowed 
to radiate, there would be a very difficult correction to the entire system of equations 
which is not of sufficient magnitude to be significant for the relatively low temperature 
and short duration of the plasmas for which these calculations are opplicable. As an 
example, the expansion could not be adiabatic if the ra<^Sation should escape. 

The seventh and last lerm is This term represents the energy in the self- 

excited oscillations of the electron cloud relative to the ions in the plasma. This is a 
difficult term to evaluate so a partial correction is employed at this time. The correc* 
tion is not large and additional corrections may be required for other plasmas. 

Limitation on this Equation of State 

The equation of state that has been presented is quite accurate until the density 
approaches solid state. The nature of the deviation is shown hy the curves in Figure 2. 
These curves show the electron density as a function of the density in grams per cm*^ 
for a constant energy added to the gas. For a perfect gas, the density of ionization of 
hydrogen gos moy be calculated from the perfect gas law with the aid of Saha*s equa- 
tion when the added energy is known. In the cases to be considered at this time, the 
added energy is either 2 eV per atom, or 5 eV per atom. The lower curve in each of 
the two graphs in Figure 2 shows the electron density as a function of the total density 
for 2 eV and 5 eV. according to the label on the graph. As the density increases to- 
ward the left in the graphs, the density of electrons continues to increase with the 
mass density . 

As a consequence of the lowering of the ionization potential by the increased 
effect of ^he microfields, the density of electrons increase faster with the mass density 
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than predicted by the perfect gas law in combination with Saha's equation. This Is 
illustrated by the curves in Figure 2. The lower curve Is the perfect gas law com- 
bined with the Saha equation. The middle curve Is for the lowering of the Ionization 
potential that is predicted by the Ecker and Kroll relotion. In every case, the gas 
Is hydrogen. 

For hydrogen gas, alone, there is a third analytical method to check the magni- 
tude of the correction and this method is precise, except the solution depends on the 
amount of energy in the other terms In the equation of state. The density of ionization 
may be calculated with the Schroedinger equation by employing the Yukawa potential. 

The Yukawa potential is for a central charge that is partially shielded, In this ca;e, by 
electrons. The calculated check on the ionization is the upper curve In the tw^ graphs 
in Figure 2. From this check calculation, it would appear that the Ecker and Kroll 
relation does not predict enough ionization. Another sink is known for the energy and 
it has been Investigated. 

The agreement between the top two curves is only fair and the agreement should be 
better. This conclusion follows from the evidence that the Schroedinger equation with 
the Yukawa potential gives an electron density that is systematically too high. The 
source of the error appears to be explainable. The Interaction energy, exists in 

the plosma to which the Schroedinger calculation is applied, but no correction for it has 
been included in the data that is plotted. When the correction for the two body inter- 
action Is included, it gives a lower electron density than the Ecker and Kroll result. At 
the same time, the curve with only the Ecker and Kroll relation should be slightly lower 
because of the three body, cluster integral correction on the energy content. The three 
body cluster integral Is mentioned below. It Is to be emphasized that the actual correction 
and the added corrections have .lOt been calculated although only arithmetic operations 
are necessary. 

In the calculation of the Interaction was found on the basis of cluster integral 

calculations on the so-called two body basis. An investlgotion showed that the Inclusion 
of three body forces woula reduce It was also shown that four body interaction would 

be exceedingly low and would have a negligible effect. In a M.S. thesis for the Mathe- 
matics Department at Oklahoma State University, Mr. H, A. Reeder calculated the three 



Figure 2. Equation of state plotted for hydrogen as a log-log relation between the 
electron der^sity vs the mass density for a total average energy of 2 eV 
per atom in the upper graph and 5 eV per atom In the lower graph. At 
the upper left in both graphs, the curve splits Into three parts to repre- 
sent the effect of lowering the ionization potential of hydrogen (a) the 
upper branch represents hydrogen as calculated by the Schroedinger 
equation with the Yukawa potential, (b) the intermediate branch by the 
Ecker and KroH relation and (c) the lower bronch contains r>o correction 
and comes from the perfect gas law and the Saha equation. 
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body cluster integral. The three body cluster integrais were for hydrogen and Reeder 
omitted the terms in the three body cluster integrals that were rK>t symmetrical. For 
ionized hydrogen, these unsymmetrical terms could not exist. His results are given 
for on added energy of 5 eV. There is o graph from his work which is designoted os 
Figure 3 in this report. The upper dotted curve is Bruce's two body calculations, the 
dohsed, intermediate curve is Reeder's cclculotion with both the two body and the 
three body terms. The lower curve is the calcuiotion from the perfect gas law and 
Soho's equation. 

PRESENTATION OF RESULTS 

With the complex equation of state that has been derived, the expansion of a 
sphere of plosmo is very different from an expansion with the perfect gas law as the 
eqijation of state. The principle effect from the extended equation of state is a conse- 
quence of term for lowering the ionization potential and the term for the inteitKil 
energy in the microfields, Before presenting the graphs, a few general comments 

are mod^ on the overall object of the computation for the expansion and on the general 
assumptions for the computation. 

General Statement of the Problem 

The cori4puter calculations ore to obtain the solution to a "gendanken" experiment. 
In this experiment, a sphere of solid aluminum is considered. At time, t = 0, each 
atom in the sphere receives the same amount of energy. In this report, there are three 
computer solutions with the odded energy per atom at 44.45 eV, 20.8 eV and 7.136 eV, 
respectively. The solution is calculated in terms of dimensionless variables. In the 
presented graphs from the solution, the dimensionless variables hove been converted to 
dimensions. The initial radius of the sphere is about 4 x 10”'^ cm. The abscissa of the 
graphs is dimensioned m terms of this starting radius ond the graphs show the growth of the 
radius as the plasma expands into a vacuum. 

In the preceding discussion, the conditions for the solution are discussed In con- 
siderable detail. One very bosic condition has been assumed without a specific state- 
ment of the postulate up to this time. This i:onditIon requires that equilibrium must exist 
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Figure 3. A log-log plot of electron density against ^he moss density for a hydrogen 
plasma. As in Figure 2, the curve divides into three branches in the jpper 
left comer of the graph. These three branches show the effect of Increas- 
ing the accuracy of the equation of state by means of a more precise calcu- 
lation of in the equation of state, (a) The upper branch shows the 
effect of the two body cluster integrals on the density of electrons, (b) The 
Intermediate bronch of the three branches shows the effect of two body cluster 
integrals plus three body, symmetrical components only, on the dertf I ty of 
electrons. This Is the most occurote curve In the group, (c) The lower 
branch in the upper left Is for the perfect gas law and the Saha equation with 
no correction for the interaction fields. 
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at ail times. In particular, this requires that ions and excited states exchange energy 
instantly. This requirement cannot possibly be correct. The mean life of an Al^^ ion 
is about 12 X 10”^ seconds. This volue was read from a curve that was part of the Ed.D 
thesis of Dr. Vernon D. Brown. Brown is presently an Associate Professor at Memphis 
State University. The mean life of aluminum ions with a higher state of ionization, 
such as Al Al^^ and Al^^" , is much shorter. This conclusion of a shorter mean life 
comes directly from the application of the uncertainty principle. The first expansion 
with an Initial energy of 44.45 eV per o^om is followed for a total time of 5.861 x 10”^ 
seconds. Ions of Al^^ and Al^^ would certainly radiate In this time interval but not 
Al^\ A correction for the mean life will certainly change the form of the solution, 
but the correction is very Involved. It Is desired, therefore, to have the solution for 
Instantaneous equilibrium. 

In the solutions that are presented, the energies are rather large. The reason for 
selecting these energies is to avoid, as much os feasible, an expansion in a region with 
densities near to solid state density. The multi -corrected equation of state doos not have 
a good correction for that region. A better equation of state Is available from the work of 
Larry J. Peery. Peery obtained the correction in his Ph.D. thesis while working with 
this group. Peery is presently a Professor at Central Methodist College at Fayette, Mo. 
His improvement in the equation of state consisted of a careful interpolation from (a) the 
accurate equation at moss densities for 10“ ‘ to 10“ ^ grams per cm'^ to (b) the density 
of solid aluminum. 

Expansion of a Plasma with an Initial Energy of 44.5 eV per Atom 

The first graph of the expansion is for an instant of time that Is very soon after the 
expansion starts. As shown in Figure 4, the exact time is 0.6845 x 10“^ seconds after 
the start of the expansion. The density is constant at the origin and is the Initial density; 
i.e., the density of solid aluminum. The average ionization potential decreases in two 
steps with a break In the slope between the two rates. The break Is near an average 
ionizotion potential of about 3 eV. There Is a corresponding small hump In the density 
at the same radius. Thic hump Is slightly masked in the reproduction. The discontinuity 
at about 3 eV indicates that this Is probably on ionization potential, or perhaps an 
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Figure 4. Density and Average Ionization Variation 
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excited state In the aluminum plasma* If this explanation is correct, more breaks 
should be observed in the graphs for other time intervals. 

The second graph in Figure 5 shows the expansion at a later time, 1 .372 x 10*^ 
seconds. The peak density has decreased by almost a factor of 10 from the density in 
the preceding graph. The average ionization potential shows two peaks, there is one 
at about 3 eV, which appears to be the same one as in the preceding g.oph. Figure 4. 
There is a second break in the slope of the average ionization potential at a little over 
1 eV. These breaks indicate that the conditions at the two radii are such as to empha- 
size an ionization potential, or an excitation level. Recall that the ionization poten- 
tials are finite and are separated. The Boltzmann relation is exponential so the curve 
is not symmetrical between the levels. 

As time continues, the expansion continues but assumes a different form. The ex- 
pansion at 2.059 x 10”^ seconds is shown in Figure 6. The density has dropped almost 
another factor of 10 from the preceding figure. The most startling difference is that the 
density decreases in the center of the sphere and the density is starting to form a cold 
shell of plasma outside the hot core. The hot and cold plasmas are identified by the 
average ionization. If the average ionization is of the order of 4.5 eV, the plasma must 
have a high Internal energy. Conversely, if the average ionization potential is nearly 
zero, there is not sufficient internal energy to be hot. There are humps in the average 
ionization which correspond to average ionization, or excitation potentials such as were 
identified in the preceding graphs. 

As the expansion continues, the formation of a hot, low density core becomes more 
clearly established. This hot core is surrounded by a cold, high density shell. This is 
demonstrated by the curves in Figure 7 which is for a time of 2.746 x 10“^ seconds after 
the expansion starts. A closer inspection of the curves in this figure show that the average 
ionization at the center of the sphere has Increased over the value at the center of the 
sphere in the preceding graph, Figu'^e 6. This Indicates that a hot core forms, and in 
additon the average ionization, or the effective temperature of the core also increases. 

In the last of the expansion curves for an initial energy of 44.45 eV, the hot core 
and the cold shell have fully formed. This expansion, at 5.861 x 10' ' seconds after 
the start is presented in Figure 8. The average ionization In the hot core is even higher 
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than the Inittol average ionization. The Initial average Ionization from the computer 
print-out Is 5.011 eV. In the hot core In Figure 8, the average Ionization is 5.550 eV* 
In this and In the preceding graphs, the average Ionization has been used as an Indlco* 
tion of the temperature. In this graph, the temperature is plotted. The temperature Is 
the curve that plunges off the bottom of the page in the center of the plot. All of the 
ordinates in this graph are linear and all of the curves except the temperature have zero 
at the bottom of the graph. In order to determine the extent of the decrease in tempera- 
ture, a plot of the temperature and the pressure are shown in Figure 9. This data is for 
the same instant of time cs Figure 8. The ordinates for both the temperature and the 
pressure are logrithmic. It is to be observed that the temperature has decreased over 
three decades at a radius that corresponds to the peak of the density in the cold shell. 

A short summary is presented of the most obvious facts that are disclosed by the 
graphs for the expansion of the very high energy plasma. 

a. The expansion forms a hot core that is surrounded by a cold shell. 

b. Humps cn the curves for the average ionization indicate the existence of the 

levels for ionization, and/or excited states in the plasma. This information was given 
to vhe con.puter by the term in the equation of state. 

c. The cold shell is really rather cold according to the computer print-out . 

The less obvious facts will be presented after comments are made on the remainder 
of the curves. There is one subject that should be mentioned. The hot core with the cold 
shell is an intrincically unstable condition. By light emission during the formation of the 
cold shell, its formation has been confirmed. Before the shell has the density range in 
Figure 8, the shell becomes unstable and breaks up to reveal the hot core. The break-up 
has been confirmed by the intensify of the emitted light. 
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Expansion of a Plasma with an Initial Energy of 20.8 eV per Atom 

For comparison with the description of the expansion of a sphere of plasma with 
very high energy content, some graphs are presented to show the expansion of a plasma 
with less than half of the high energy in the first expansion. The actual energy con- 
tents are 20.8 eV per atom In contrast to 44.45 eV per atom. The same number of atoms 
are involved in each case. 

The first graph for the expansion is given for a time of 1 . 108 x 10’^ seconds aftc r 
the start of the expansion. The curves in this graph, Figure 10, are for the lower energy, 
20.8 eV per atom. They are to be compared with the high energy graphs that are pre- 
sented In Figure 4 and Figure 5 which are for 0.6845 x 10”^ seconds and 1 .372 x 10"^ 
seconds, respectively, after the expansion started. The curves for the lower temperature, 
in Figure 10, show a slower expansion of the high density material in the core. This 
would be expected. The breaks in the curves for the averoge ionization ere more pro- 
nounced than at the higher temperatures. A break in the average ionization still occurs 
near the average value of 3 eV. The exact value fluctuates about the 3 eV level. This 
is to be anticipated and depends on the density, the temperature, the location of excita- 
tion levels and other factors. 

As the expansion continues, the density still corresponds to the initial solid state 
density although the expansion started 2.261 x lO""^ seconds before the curves in Figure 11 . 
The '*knees" in the average Ionization curve are more pronounced than ever. 

As the expansion continues, the start of the formation of the hot, low density core 
is well defined in Figure 12. This graph illustrates the state of the expansion and of the 
average ionization after 3.413 x 10“‘* seconds. There is nothing particularly new with 
respect to the results in other graphs. 

The hot, low density co/e is well defined in the next graph. The set of curves in 
Figure 13 were obtained 4.565 x 10*^ seconds after the start of the expansion. The curve 
for the average ionization looks rather distorted from the presence of strongly emphasized 
ionization and/or excitation levels near 3.0 eV and below 1 .0 eV. This is approximately 
the relative distribution between hot core and cold shell at which the combination be- 
comes so unstable in the practical case that It breaks up, or more probably, It just de- 
composes. This is a highly arbitrary conclusion, although there are a few threads of 
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evidence from the light emission. These are the experiments that have been reported 
by Brown. 

The last graph for the expansion of the plasma with an Initial energy of 20.8 eV 
per atom is designated as Figure 14. This graph shows the expansion at 9.808 x 10“^ 
seconds after the start of the expansion. An observation of the graphs Indicate that the 
expansion has been computed far beyond the limit at which a practical plasma would be 
expected to become so unstable that it would disintegrate. The average Ionization curve 
has become rather smooth but It still shows evidence of the 3.0 eV and the 1 .0 eV Ioni- 
zation and/or excitation levels. 

In comparing the graphs for an Initial energy of 20.8 eV per atom with the preceding 
graphs for an Initial energy of 44.45 eV per atom^ It Is difficult to select the most Impor- 
tant factor on which to concentrate. A very interesting subject appears to be the approxi- 
mate mean radius of the expanding sphere when the density in the core has just dropped 
into the range from 0.1 to 0.2 x 10“" grams per cm^ . For the higher Initial energy, this 
radius is in the range between Figure 7 and Figure 8. The subjectively chosen value is 
0.0027 cm. In a similarly unreliable manner, the radius is chosen just later than the in- 
stant that is represented In Figure 13. Then for the initial energy of 20.8 eV per atom, 
the value of the radius is taken to be 0.00202 centimeters. Without reference to the time 
that is required to expand to these limits. It is observed that roughly similar conditions 
exist between the hot core and the cold shell at the time that these radii were selected. 
These values are selected while the curves are in close proximity to this discussion In 
order that they may be easily inspected. There will be a discussion of these measurements 
after the next set of graphs is presented. 

Expansion of a Plasma with an Initial Energy of 7. 136 eV per Atom 

In order to acquire more Information on the effect of the total Initial energy on the 
expansion of a sphere of plasma, some curves are shown for the case when ♦he total Initial 
energy is 7.136 eV per atom. This energy is less than 1/6 of the total energy for the case 
when the Initial energy was 44.45 eV per atom. 

With this lower energy, the most Important question is to determine If the cold shell 
still forms around a hot core and the approximate time that is required for this shell to be- 
come roughly formed. The curves in Figure 15 represent the state of formation at 4.636 x 10 
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seconds after the start of the expansion when the initial energy was 7. 136 eV per atom. 

At this time, the hot core is starting to form and the dense, cold shell is beginning to 
narrow into a peak. The density, as reported by the scale of the ordinate, is down by 
a factor of about 20 from the initial density. 

As the exponsion continues, the curves in Figure 16 show the expansion at 7.059 x 
lO""^ seconds after the start of the expansion. The initial energy of the plasma was 7.136 eV 
per atom. The curve for the average ionization clarifies some of the up and down variation 
in the position of the ionization and/or excitation levels near 1.0 eV. In other curves, the 
level sometimes appeared to be above 1 eV and sometimes below that level. Actually, 
there are three levels. One is near 1 .3 eV, a minor level Is near 0.9 eV and a very strong 
level is at 0.6 eV. The peak of the density in the cold shell is desired for the same pur- 
pose. 

DISCUSSION OF THE SIGNIFICANCE OF RESULTS FROM THE EXPANSION OF PLASMAS 


Some superficial results are repeated to start this discussion. These superficiol results 
were collected at the end of the presentation of the curves on the expansion of the plasma 
with an initial energy of 44.45 eV per atom. 

a. The expansion forms a hot core that is surrounded by a cold shell . 

b. Humps, or breaks, in the curvature of the curves for the average ionization indi- 

cate the presence of levels for ionization and/or excitation in the plasma. This informa- 
tion was given to the computer by the term that is designated as corrections to 

the perfect gas law. These le^'els are not fixed but may appear to be raised, or lowered by 
the effect of other terms in the equation of state. 

c. The curve for the temperature in Figure 9 shows that the :>o-called cold «nell is 
really cold. In this figure, the temperature falls over three decades within a radial dis- 
tance that is less than the width of the cold shell . 

After the preceding list of specific comments, there are two general subjects which 
require some discussion. The first subject deals with the formation of the cold shell. It 
involves the means by which the high temperature and energy in the initial sphere is con- 
verted to a high radial velocity. After the comments on the cold shell, the second problem 
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for discussion is the probable evidence that the hot core is squeezed and slowed in 
expanding. 

In the consideration of the first subject on the conversion of temperature to kinetic 
energy ; the following facts are emphasized. There is no energy added to, or subtracted 
from the sphere during its "free" expansion into a vacuum. There is only a change in the 
nature of the energy storage. Although the expansion into the vacuum is "free", there 
is a change in entropy as decreases to zero with the expansion. It is to be recalled 
that is the energy in the microfields as calculated by cluster integrals. At the start 
of the expansion, all of the sinks of energy, v<hich are indicated in the equation of state, 
are filled and the values are placed in the computer which is run until the distribution con- 
verges on the correct values. In this way, the energy distribution is in equilibrium. When 
expansion is ol lowed to start, the ions and electrons on the outside of the sphere start to 
move outward from the sphere, but only a small part go out radially. As time progresses, 
the motion of the outermost ions and electrons become more and more directed outward along 
radial lines. At the start of the expansion, all of the kinetic energy is randomly directed. 
With randomly directed velocities, the temperature is given by the relation. 


f '"ion I 

where m^ is the mass of an electron, and is the mean square velocity. and 

are the same values for the ions and k is the Boltzmann constant. 

By ♦he formulation of the equation of state, the radially directed components of the 
plcsrna flow i> no longer recorded as a temperature. The temperature is composed only of 
the randomly directed velocities. As an example, the rate of expansion approaches a 
consr.7nr as the sphere grows larger. This constant velocity is obtained when all of the 
energy in the surface atoms is directed radially outward. As a consequence of the dis- 
appearance of the random velocities, the temperature falls to very low values as is appar- 
ent by the c^rve for the temperature in Figure 9. 

The second subject for consideration is the relative rotes of expansion in the early 
stages while the hot core Is being formed. The radius of the expanding sphere was esti- 
mated for the time that the hot core and cold shell were just barely formed. At the end 
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of the presentation of the graphs for on initial energy of 20.8 eV, there is a collection 
of subjective data. This collection of data is to show the radius of the peak of the 
dense, cold shell at a particular instant in the expansion. It is the instant at which the 
hot core and cold shell are just developed and the lowest density of the hot core has just 
reached the range between 0.1 and 0.2 x 10”^ grams per cm® . For an initial energy of 
7.136 eV per atom, the radius is taken to be 1 .7x 10“® centimeters. This data on radii 
is fragmentary and very subjective. There has been hesitation to add it to this report. 
The data is presented in Table III. 


TABLE III 


Subjective Data Illustrating Self-Pressure of Hot Core 


Initial 

Energy 

Tlnitial 
^ Energy 

Ratio 
llnitial 
\ Energy 
■ 2.67 

r 

cm X 10~® 

Predicted r 
r X Ratio 

44.45 

6.66 

2.50 

2.7 

6.75 

20.8 

4.56 

1.71 

2.02 

3.45 

7.136 

2.67 

1. 

1.7 

1.7 


In Table III, the columns have the following significance. The first column is ob- 
vious. The second column is the square root of the first column which means that the 
terms in this column are roughly proportional to the velocity of the final rate of ex- 
pansion. The third column ranks the relative rates of expansion. They are ranked 
relative to each other by means of the lowest energy term; i.e. relative to an initial 
energy of 7.136 eV per atom. The fourth column gives ihe subjectively selected radii 
with no consideration of the instant of time after the start of the expansion at which the 
radius was attained. The fifth column is the product of the subjectively selected value 
of r by a quantity that is proportional to the final rate of expansion. 

The values in the right hand column are the radii that would be expected if no self- 
pressure was exerted on the hot core by the lowering of the ionization potential and by 
the energy in the term This data is interesting but it is certainly not con- 

clusive. In trying to interpret these results, it should be considered that the velocity of 
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expansion starts slowly and Increases with time to the predicted final rate. 

There are two more subjects which should be mentioned briefly. The average 
ionization at the center of the sphere is always higher in the last, or the last few 
graphs, than the average ionization at the start of the expansion. This is rather 
surprising in consideration of the lower density and the corresponding lower pressure. 

No further discussion of this subject is given in this report. There is another interest- 
ing fact. The pressure falls very precipitously just inside the cold shell. This pressure 
is not a proper scaler pressure. It is evidence that a tensor should represent the pressure. 
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PHASE THREE 

ADAPTATION OF A COMPUTER PROGRAM TO A SOLUTION FOR METEOROID 

BUMPER PENETRATION 

INTRODUCTION 

The objective of this study is to develop a reliable computer program for the 
penetration of a thin plate of aluminum by a sphere of quartz. The thin plate of 
aluminum simulates a bumper plate for micrometeoroids. A good equation of state is 
available ^or quartz as well as for aluminum. As the aluminum equation of state was 
derived, the equation for quartz is based on an experimentally determined Hugoniot 
curve. Then an equation of state of the Mie-Grunesen type was derived. In addi- 
tion, an adjustment in the computer program permits the quartz to have any desired 
porosi ty . 

At the start of this study, a computer program was available for the formation of 
a crater by the impact of a porous sphere of rock at hypervelocity onto a thick slab of 
aluminum. Mathematically, the thick slab is a semi-infinite slab of aluminum. The 
available computer program was formulated by Dr, B. A. Hardage, who is an ex-member 
of this group. Dr, Hardage is presently a member of the Geophysical Group for Phillips 
Petroleum Company. Hardage employed this program with considerable success to de- 
termine the size and shape of craters that were formed by hypervelocity impact. As will 
be indicated, it was suitable for that purpose, but simplifying assumptions by Hardage 
were not appropriate for the penetration of a thin plate of aluminum. With his program, 
incorrect results were found for the penetration of a thin plate. 

The modifications to Hardage 's program are being made entirely by Mr. Mark J. 
Hooker. The modifications have proved to be much more extensive than was anticipated. 
As a consequence, the final computer progrom is considered to be developing Mr. Hooker's 
capabilities far beyond those of the average student who writes a thesis as part of the 
requirements for an M.S. degree. Since this is Hooker's first major computer program, 
he has had considerable assistance beyond that which is received by the usual student. 

Dr, R. E. Bruce helped in changing the differencing method. Dr. Hardage has given 
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major assistance in emphasizing the limitations of his proyram, disclosing the exact 
nature of some assumptions that were nor clear, and approving some of the modifica- 
tions. 

The following discussion presents the basis for the changes in the computer pro- 
gram and comments on the results from these changes. Since the modifications vary so 
widely in importance and in the difficulty of the actual change, they will not be listed 
separately in this introduction. The changes v/ill be grouped under headings which will 
indicate the basis for some of the corrections. 

TECHNIQUE TO DETERMINE MOMENTUM AT BOTTOM SURFACE OF THIN PLATE 

For general information and to assist in formulating the remainder of the problem, 
it is desired to be able to determine the momentum that reaches the bottom of the thin 
plate. The only part of immediate Interest is from the initial shock when the sphere of 
rock is incident on the thin plate. The method of obtaining this information Is described 
to assist in interpreting the graphs, ii hcs not been changed for it serves another p 
pose, it decreases the cost of trouble shooting other changes in the overall program. 

This last use is discussed before the momentum is considered. 

In the curves that have been presented and in the one curve that is given in this 
report, no tension force is assumed to be exerted between the curvilinear cubes with 
which the computer deals. Compression *vas always and still is in the computer program. 
This is a hold-over from the Hardage program. His program dealt only with compression 
in the formation of a crater. The spray that "flew up " from the impact zone was re- 
leased if the internal energy exceeded the energy to melt aluminum. These instructions 
have not been changed up to the present time. In the portion of the problem in which 
tension predominates, it permits the plate to disintegrate. There have been so many 
problems in the compression portion of the problem that this characteristic has been an 
aid. One advantage of this omission of rupture forces between the curvilinear cubes is 
the considerable reduction In computer time and in computer printout while correcting 
other features In ihe computer program. 

A second, more Important advantage is that this technique, with p-^oper sc re- 
guards, may be employed to estimate the momentum that Is involved in the reflection 
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of the initial shock from the back face which is the bottom face J he thin plate. 

This may be illustrated by reference to the graph In Figure 17. On the bottom side of 
the thin plate, there is a hemispherical gap of missing material between columns 10 
and 14, counted from right to left. Recall that the Initial bottom face was horizontal. 

This region shows a part of the effect of the incident momentum of rr.e sphere of rock. 

When the rock reaches the plate, a shock starts from the initial region of contact between 
the sphere and the plate. The portion of this shock front, which enters the plate, reaches 
the rear face and “snaps off“ a chunk of aluminum as It is reflected from the rear face. 

The missing material is the mass which is compressed as the Initial shock passes it on the 
path to the back face. This shock compression gives the material that energy which is 
required to exceed the energy to melt the material. The velocity of departure of this 
rock may be measured by observation of the graphs from cycles which precede the cycle 
that is presented in Figure 17. Witli the mass and the velocity known, the momentum is 
determined. This is an application of a simple, familiar experiment in physics. Suspend 
ten balls in a line and in contact. Lift one ball and let it impact on the second ball . 

Only the tenth ball, at the end of the line, flies up. 

ENERGY FOR HYPERVELOCITY SEPARATION OF FRAGMENTS 
OF THIN ALUMINUM PLATE 

In the preceding section, the force to rupture the aluminum was mentioned. It was 
also indicoted that no value for the strength of aluminum had been inserted in the pro- 
gram up to this time. The computer program must always consider that the break occurs 
over the surface of a curvilinear cube. The failure may be rupture in tension, or it may 
be a failure In shear. There Is some evidence in small leglons thot the failure me/ approxi- 
mate to a plane; however, the most of the ruptures are very rough. This appears to elimi- 
nate shear as c major factor. From the observations, most of the ruptures appear to be in 
failures in tension. The strength in tension for rapid breaks is difficult to obtain, or 
recognize among the available data. The coruputer calculation for the graph in Figure 17 
was based on the assumption of an Impact at 7.2 kilometers per second. 

For practical purposes, a simple assumption 'v;li be made in order to obtain a good 
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Figure 17. A Topical Graph of the Performance of the Computer Program for the Impact 
of a Sphere ^ock ^Quartz^ on a Thin Plate o*^ Aluminum. The figure rep* 
resents -e si c^'-er obout 2O0 c>xles or the computer. Hold the graph 

with The ''or-o.ves' squares at the top. The r^^ht bond side of : tC graph rep- 
resents the position of the diometer of the sphere at the start of the Impact. 
The sphe r me es *rom. the top ‘owerd the bottom o*' 'he graph. At the start 
of the im::c:c t'^e hemisphere extends from 1 to 20 with a radius of 10 spaces 
in the ^errical direction. The type and the amount of the information is indi- 
cated b/ the following description of the symbols. Only about half of the 
information per curvilinear cube is presented Omitted information includes 
density, internal energy, kinetic energy, radial momentum, transverse 
momentum, etc. 

In this graph, the symbols are listed below. The symbol in each curva- 
linear cube. 

o Partly filled with rock, remainder vacuum 

O Completely filled with rock 

X Partly filled with aluminum, remainder vacuum 

X Cdmpletely filled with aluminum 

II Portly filled with aluminum, partly filled with rock 

III Portly filled wth aluminum, partly filled with rock and 

partly filled with vacuum 

The fractional volume of aluminum and of rock is given in eoch case. 




or 
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set of values for the final solution. There are two variables to be considered, the rate 
of separation of the parts and the temperature. High rates of separation will reduce the 
amount and perhops even the possibility of ductile failures by the flow of the metal, or 
liquid. The brittle failures are not of the type which occurs when a crack enters a highly 
stressed region so the crack constitutes a stress “riser'*. It is more likely that the pi^ $ 
will be simply jerked, or knocked apart. Without considering the effect of the rate in 
more detail, the effect of temperature is considered. 

From the internal energy, that is in the print-out from the computer, the tempera- 
ture of the metal, or liquid will be known. Below the melting point, which is only 659. 7^C, 
the strength of aluminum decreases in a manner that is readily available from references. 

These values will probably be increased in the final solution provided the effect of the high 
rate of separation is introduced. For the first, approximate values, the liquid will be re- 
quired to withstand some force in tension. T*^e exact value will be selected at the last 
minute. Above the boiling temperature of aluminum at 2057^C, the aluminum is a plasma 
and will “fly“ apart. 
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APPENDIX A 

Typical Substitution of Experimental Results to Evaluate Constants 

The evaluation of the constants by employing experimental data illustrates some features 
of the equation. Por this purpose, an evaluation is given below. The general equation, 
Equation 10, has the following form 



(Pa)° 9P,ro -7P, 
24 (PaVp.ro)^ 



The known constants and the known data from one experiment may be inserted to partially 
evaluate the constants, Pj, Pg and P,. 

Pi = 7.5313 X 10® D 

s^ = 8 cm 

The dependence of the second constant, Pg on various values is given by the relation 
Pg=4p - 3XPi 

i is assumed for this evaluation to have the value, 0.01 
0 is the density of aluminum, 2.675 

Oj is the average density of the stacked sheets of paper, 0.4609 
Then the numerical value of Pg is given by the relation 
Pg =11.060 - 0.01383 = 11.0462 

The dependence of the third constant, P..,, on various values is given by the relation 


Pa = -6 i S 

Z is cisumed for this evaluation to have the value, 0.01 
S is taken to be 10 grams per cm 
Then the numerical value of Pa is given by the relation 


Pa = -0.6 
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Substituto these constants into the equation at the start of this Appendix. 


O T co,o _(Pa)M9x 11.0462 x 0.09525 + 4.2) , 7 /P,f 

8 X 7.5313 X 10 D - ^ " 24^^) 


5.025 X I0"'D = 


3.559:5 


+ Pa>-0? 24 ■ 24 \P 

Continue with the evaluation by substituting for the remaining values of and P. 


. AOC ,0.0 (11.0462)® 13.66935 ^ 

6.025x 10* D 6 + 11.0426 x 0.09525)'^ 2? 


1 _ ( 

24 \ -0. 


0462 r 


Two features of this numerical relation should be noted. The first is the critical nature of the 
denominator in the first term. The denominator is the difference between two relatively small 
numbers. The difference is positive for this set of values. The other feature to note is that, 
in the second term, the numerical value inside the bracket is negative. Since this negative 
value is raised to an even power, 8, the value of the second term is always positive. 

z AOC (11.0462)® 13.66535 ^ 7 /11.0462')® 

6.025 X 10* D - - + ^ [-ITT-) 

Continue the simplification of this relation 

z AAC ,A10^ _ (11.0462)® 13.66535 ^ 7 /'ll. 0462)® 

X It u (0.45215)'^ 24 24 V 0.6 / 

After a few arithmetical steps, the equation reduces to 


D =2.716 



